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Abstract 

Using a loop-cluster algorithm we investigate the spin | Heisenberg 
antiferromagnet on a square lattice with exchange coupling J and an 
additional four-spin interaction of strength Q. We confirm the existence 
of a phase transition separating antiferromagnetism at J/Q > Jc/Q 
from a valence bond solid (VBS) state at J/Q < Jc/Q- Although our 
Monte Carlo data are consistent with those of previous studies, we do not 
confirm the existence of a deconfined quantum critical point. Instead, 
using a fiowgram method on lattices as large as 80^, we find evidence for a 
weak first order phase transition. We also present a detailed study of the 
antiferromagnetic phase. For J/Q > Jc/Q the staggered magnetization, 
the spin stiffness, and the spinwave velocity of the antiferromagnet are 
determined by fitting Monte Carlo data to analytic results from the 
systematic low-energy effective field theory for magnons. Finally, we 
also investigate the physics of the VBS state at J/Q < Jc/Q, and we 
show that long but finite antiferromagnetic correlations are still present. 
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1 Introduction 



Undoped antiferromagnets, which can be modeled with the spin | Heisenberg Hamil- 
tonian, are among the quantitatively best understood condensed matter systems. To 
a large extent this is due to an interplay of the very efficient loop-cluster algorithm 
[H, 0,l3| with the effective field theory for antiferromagnetic magnons 0, [g, 0, @] • 
In particular, applying chiral perturbation theory — the systematic low-energy effec- 
tive field theory foi^Goldstone bosons — to antiferromagnetic magnons, Hasenfratz 
and Niedermayer |9| have derived analytic expressions for the staggered and uni- 
form susceptibilities. By comparing these expressions with very accurate Monte 
Carlo data obtained with a loop-cluster algorithm, the staggered magnetization, 
the spin stiffness, as well as the spinwave velocity of the Heisenberg model have 
been determined very precisely 0]. In particular, the resulting values of these 
low-energy parameters are in quantitative agreement with experimental results on 



undoped antiferromagnets [10 



High-temperature superconductors result from doping their antiferromagnetic 
precursor insulators. With increased doping, antiferromagnetism is destroyed be- 
fore high-temperature superconductivity emerges. Understanding the doped systems 
from ffist principles is very difficult because numerical simulations of microscopic 
systems such as the Hubbard or t-J model suffer from a severe fermion sign problem 
at a non-zero density of charge carriers. In the cuprates, antiferromagnetism and 
high-temperature superconductivity are separated by a pseudo-gap regime. It has 
been conjectured that this regime is connected to a quantum critical point with un- 
usual properties. In particular, the Neel order of the antiferromagnet may give way 
to a spin liquid phase without long-range magnetic order before the phase coher- 
ence of the Cooper pairs of high-temperature superconductivity sets in at somewhat 
larger doping. 

According to the Ginzburg-Landau- Wilson paradigm, a direct phase transition 
separating one type of order from another should generically be of ffist order. This 
paradigm has recently been challenged by the idea of deconfined quantum criticality 
[III . A deconfined quantum critical point is a second order phase transition directly 
separating two competing ordered phases as, for example, an antiferromagnet or a 
superfluid from a valence bond solid (VBS). There are two types for VBS order: 
columnar and plaquette order, which are illustrated in figure 1. At a deconfined 
quantum critical point, spinous — i.e. neutral spin ^ excitations which are confined 
in the two ordered phases — are liberated and exist as deconfined physical degrees 
of freedom. It was conjectured that the continuum field theory that describes a 
deconfined quantum critical point separating an antiferromagnet from a VBS state 
should be a (2 -|- l)-dimensional CP(1) model with a dynamical non-compact U{1) 
gauge field. This theory is expected to be in the same universality class as an 
0(3) non- linear a-model in which the creation or annihilation of baby-Skyrmions 



is forbidden 12| . The resulting conserved number of baby-Skyrmions gives rise to 
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Figure 1: Columnar (a) and plaquette (b) type of VBS order. The solid bonds 
indicate groups of spins that preferentially form singlets. 



an additional U{1) symmetry. In ISj it has been argued that — upon doping — 
a deconfined quantum critical point separating antiferromagnetism from VBS order 
may extend to a spin liquid phase, thus providing a possible explanation for the 
pseudo-gap regime in under-doped cuprates. 

Establishing the existence of deconfined quantum criticality in an actual physical 
system is a non-trivial issue. For example, numerical simulations of microscopic 
models with a transition separating superfluidity from VBS order found a weak 



first order transition |l4j . Il5l |. Using their fiowgram method, in a detailed study of 
another superfluid-VBS transition Kuklov, Prokof'ev, Svistunov, and Troyer have 
again established a weak first order transition instead of a quantum critical point 



161], thus confirming the Ginzburg-Landau- Wilson paradigm also in that case. 



While unbiased numerical simulations of sufficiently strongly doped antiferro- 
magnets are prevented by a severe fermion sign problem, Sandvik has pointed out 
that there is no sign problem in the spin ^ Heisenberg model with a particular 



four-spin interaction similar to ring-exchange [17[. For this model, he presented 



numerical evidence for a deconfined quantum critical point separating antiferromag- 
netism from VBS order. The quantum Monte Carlo study of I TII w as performed 
using a projector Monte Carlo method in the valence bond basis isl, 19] and was lim- 
ited to zero temperature and to moderate volumes. Recently, Melko and Kaul have 
simulated the same system on larger lattices at finite temperature using a stochastic 



series expansion method 20|]. Both studies 17|, l20l] conclude that the transition 
belongs to a new universality class that is inconsistent with the Ginzburg-Landau- 
Wilson paradigm. As we will discuss, this conclusion rests on the use of sub-leading 
corrections to scaling, which can, however, not be determined unambiguously from 
the data. In this paper, we apply a rather efficient loop-cluster algorithm to the 
same system. This has allowed us to also reach large volumes. In order to decide 
if the transition is second or weakly first order, we have implemented the fiowgram 
method of 16||. Our data provide evidence that the transition is weakly first order. 
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i.e. the Ginzburg-Landau- Wilson paradigm is again confirmed. This means that an 
SU{2) or U{1) invariant system, for which the phenomenon of deconfined quantum 
criticality can be firmly established, has yet to be found. Finding such a system 
is non-trivial, in particular, since it should be accessible to accurate first principles 



numerical simulations. In this context, it is interesting to consult [2l|, |22 



Besides studying the phase transition, we also investigate in detail how antiferro- 
magnetism is weakened. In particular, we extend the results of 0, 0] by determining 
the staggered magnetization, the spin stiffness, as well as the spinwave velocity of the 
Heisenberg antiferromagnet as functions of the strength of the four-spin interaction. 
In addition, we investigate some properties of the VBS phase. 

The rest of the paper is organized as follows. The Heisenberg model with four- 
spin interactions as well as some relevant observables are introduced in section 2. In 
section 3 the weakening of antiferromagnetism is studied by comparing Monte Carlo 
data with analytic predictions from the systematic low-energy effective theory for 
magnons. In section 4 the phase transition and, in particular, the question of its 
order is investigated. Some properties of the VBS phase are studied in section 5. 
Finally, section 6 contains our conclusions. 



2 Spin Model and Observables 

In this section we introduce the microscopic Heisenberg Hamiltonian with four-spin 
interaction, as well as some relevant observables. 



2.1 Heisenberg Model with Four-Spin Interaction 

Let us consider the spin | Heisenberg model on a 2-dimensional periodic square lat- 
tice of side length L with an additional four-spin interaction defined by the Hamil- 
tonian 



- - 1 - - 1 

{Sx ■ S^^l — -)(>S'^+2 ■ ~ 7) 



4' 



(2.1) 



Here = \Bx is a spin | operator located at the lattice site x and i is a vector of 
length a (where a is the lattice spacing) pointing in the i-direction. The standard 
exchange coupling J > favors anti-parallel spins. The four-spin coupling Q > 



^Renormalization group arguments have been used to demonstrate the existence of deconfined 
quantum criticahty in an SU (A^)-invariant system at very large N . 
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favors the simultaneous formation of singlet pairs on opposite sides of an elementary 
plaquette. Sandvik has pointed out that quantum Monte Carlo simulations of this 



four-spin interaction do not suffer from the sign problem [171]. Indeed, the fact that 
it can be treated reliably in numerical simulations is the main reason to consider 
this particular form of the coupling. 



2.2 Observables 

Obviously, the above Hamiltonian commutes with the uniform magnetization 

M = Y,Sx. (2.2) 

The order parameter for antiferromagnetism is the staggered magnetization 

Ms = ^(-l)(^^+^'^)/"S;. (2.3) 

x 

A physical quantity of central interest is the staggered susceptibility 

Xs = j2iyt {M!{0)M^{t)) = ^iyt ^Tr[M^iO)M^{t)exp{-pH)], (2.4) 

the integrated correlation function of the 3-component of the staggered magnetiza- 
tion operator. Here (3 is the inverse temperature and 

Z = TTexp{-(3H) (2.5) 

is the partition function. Another relevant quantity is the uniform susceptibility 

Xu = j2jy^ {M\m\t)) = j2jy^ |Tr[M3(0)M3(t)exp(-/5i/)], (2.6) 

the integrated correlation function of the uniform magnetization. Both Xs and Xu 
can be measured very efficiently with the loop-cluster algorithm using improved esti- 
mators j2| • In particular, in the multi-cluster version of the algorithm the staggered 
susceptibility 

X. - (i: 

is given in terms of the cluster sizes \C\ (which have the dimension of time). Similarly, 
the uniform susceptibility 
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is given in terms of the temporal winding number Wt = X]c^<(^)' which is the 
sum of winding numbers Wt{C) of the loop-clusters C around the Euclidean time 
direction. In complete analogy, the spatial winding numbers Wi = Wi{C) define 
two spatial susceptibilities 



These susceptibilities measure the response of the system to a twist in the spatial 
boundary conditions. 

A natural order parameter that signals a VBS state is 



In a VBS state with columnar order either Di or D2 has a non-vanishing vacuum 
expectation value. In a VBS state with plaquette order, on the other hand, one of 
the linear combinations Di ± D2 has a non-zero expectation value. In numerical 
simulations, it is easier to investigate an alternative pair of order parameters which 
just count the number of spin flips in the configurations contributing to the path 
integral. We define the order parameter Di as the difference between the number of 
spin flips on nearest-neighbor bonds in the i-direction with an even and an odd value 
of Xi/a. It should be noted that such flips can be due to both the standard two-spin 
coupling of strength J and the four-spin coupling of strength Q. The corresponding 
probability distribution p{Di,D2) is useful for investigating the nature of the VBS 
state. 



3 Weakening of Antiferromagnetism 

In this section we investigate the weakening of antiferromagnetism. First, we briefly 
review some results of the systematic low-energy magnon effective field theory. Then 
Monte Carlo data obtained with a loop-cluster algorithm are used to determine the 
values of the low-energy parameters of the effective theory. 

3.1 Low-Energy Effective Theory for Magnons 

The low-energy physics of antiferromagnets is determined by the SU{2)s spin sym- 
metry which is spontaneously broken down to U{l)s. As a result, there are two mass- 
less Goldstone bosons — the antiferromagnetic spinwaves or magnons. Chakravarty, 
Halperin, and Nelson jj] were first to describe the low-energy magnon physics by an 
effective field theory — the (2 + l)-d 0(3)-invariant non-linear a-model. In analogy 




(2.9) 




(2.10) 



X 
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to chiral perturbation theory for the pseudo-Goldstone pions in QCD, a systematic 
low-energy effective field theory for magnons was developed in 0, EBB- The 
staggered magnetization of an antiferromagnet is described by a unit-vector field 

e(a;) = (ei(x), 62(0;), 63(0;)), e{xf = 1, (3.1) 

in the coset space SU{2)s/U{l)s = S"^. Here x = {xi,X2,t) denotes a point in 
space-time. To leading order, the Euclidean magnon effective action takes the form 

S[e\ = j d'x dt ^ (^d,e- d,e+ ^dte- dte^ . (3.2) 

The index i G {1,2} labels the two spatial directions, while the index t refers to 
the Euclidean time-direction. The parameter ps is the spin stiffness and c is the 
spinwave velocity. At low energies the antiferromagnet has a relativistic spectrum. 
Hence, by introducing xq = ct the action takes the manifestly Lorentz-invariant 
form 

S[e\ = J d^x dxo ^d^e- d^e. (3.3) 

The ratio ^ = c/ (27rps) defines a characteristic length scale which diverges when 
antiferromagnet ism disappears at a second order phase transition. 

Hasenfratz and Niedermayer have performed very detailed calculations of a va- 
riety of physical quantities including the next to next to leading 2-loop order of 
the systematic expansion For our study their results for finite temperature and 
finite volume effects of the staggered and uniform susceptibilities are most relevant. 
Depending on the size L of the quadratic periodic spatial volume and the inverse 
temperature (3, one distinguishes cubical space-time volumes with L (3c from cylin- 
drical ones with (3c ^ L. The aspect ratio of the space-time box is characterized 
by 

In the cubical regime the volume- and temperature-dependence of the staggered 
magnetization is given by 

. ^ |l + + (_±.)^ + 3A(/)] +...}, (3.5) 

where A^s is the staggered magnetization density. The uniform susceptibility takes 
the form 




(3.6) 



The functions /?«(/), [3i{l), and ip{l) are shape coefficients of the space-time box 
defined in 0]. 
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Q/J 


I3J 


L/a 


XsJa 




0.1 


20 


34 


743.0(1.8) 


8.285(26) 


0.1 


20 


36 


829(2) 


9.312(25) 


0.5 


16 


42 


625.1(1.6) 


7.561(20) 


0.5 


16 


44 


683(2) 


8.310(20) 


1 


12 


42 


333.7(1.3) 


5.361(20) 


1 


12 


46 


396.3(1.5) 


6.434(22) 


2 


12 


66 


470.8(1.6) 


5.598(20) 


2 


12 


68 


497.7(1.6) 


5.960(22) 


3 


10 


78 


383.8(1.2) 


5.310(21) 


3 


10 


82 


420.8(1.2) 


5.914(22) 


4 


10 


94 


415.9(1.5) 


5.011(23) 


4 


10 


96 


431.4(1.5) 


5.229(26) 



Table 1: Some numerical data for the staggered susceptibility and the temporal wind- 
ing number squared {WD obtained with the loop-cluster algorithm. 

3.2 Determination of the Low-Energy Parameters 

We have performed numerical simulations of the Heisenberg model with four-spin 
interaction for a variety of lattice sizes L/a ranging from 24 to 112 at inverse tem- 
peratures between /3J = 10 and 20. Remarkably, just hke the ordinary two-spin 
couphng, the additional four-spin coupling can also be treated with an efficient 
loop-cluster algorithm. The algorithm, presently implemented only in discrete time, 
will be described elsewhere. All simulations described in this section have been per- 
formed at three different lattice spacings in discrete time, which allows a reliable 
extrapolation to the continuum limit. Some numerical data (extrapolated to the 
time-continuum limit) are listed in table 1. For fixed J and Q all data for Xs and Xu 
have been fitted simultaneously to eqs. (13.51) and (13.61) by using the low-energy con- 
stants A^s, ps, and c as fit parameters. The fits are very good with x^/d.o.f. ranging 
from 0.5 to 2.0. Typical fits are shown in figures 2a and 2b. The corresponding 
results are summarized in table 2 and illustrated in figures 3 and 4. One observes a 
substantial weakening of antiferromagnetism. In particular, as one goes from (5 = 
to Q = 4J, the staggered magnetization M.^ decreases by a factor of about 3, while 
the correlation length = c/{2Txps) increases by a factor of about 5. Interestingly, 
in units of J, the spin stiffness ps is more or less constant. The increase of ^ with 
Q is thus due to an increase of the spinwave velocity c (in units of Ja^). When 
antiferromagnetism disappears at a second order phase transition, the correlation 
length ^ diverges. This is possible, only if ps goes to zero at the transition. Since 
the system interacts locally, any excitation travels with a finite speed, and hence c 
cannot go to infinity. In the next section we will present numerical evidence for a 
ffist order phase transition. In that case, ps remains finite at the transition. 
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Figure 2: Fit of the finite- size and finite-temperature effects of the staggered suscep- 
tibility Xs (a) o-nd the temporal winding number squared (Wf) (b) to results of the 
effective theory in the cubical regime for various values of Q/ J . 
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Q/J 




Ps/J 


c/(Ja2) 







0.3074(4) 


0.186(4) 


1.68(1) 


1.44(3) 


0.1 


0.2909(6) 


0.183(6) 


1.88(3) 


1.64(3) 


0.5 


0.2383(7) 


0.182(6) 


2.73(4) 


2.39(4) 


1 


0.1965(7) 


0.194(7) 


3.90(6) 


3.19(5) 


2 


0.149(1) 


0.194(9) 


5.98(14) 


4.91(12) 


3 


0.122(1) 


0.192(8) 


7.97(16) 


6.60(14) 


4 


0.106(1) 


0.218(13) 


10.50(31) 


7.67(26) 



Table 2: Results for the low-energy parameters M.s, Ps, and c as well as the length 
scale ^ = c/{2ttPs) obtained from fitting Xs o-nd Xu to the analytic expressions of 
eqs. 113. 5\) and /13. 6\) from the magnon effective theory. 




Figure 3: The staggered magnetization Aig a function of Q/J, obtained from the 
fits to the magnon effective theory results for Xs and Xu- 
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Q/J 



Figure 4: The length scale ^ = c/{2nps) as a function of Q/J, obtained from the fits 
to the magnon effective theory results for Xs (^nd Xu- 

4 Phase Transition between Antiferromagnetism 
and VBS Order 

In this section we study the phase transition at which antiferromagnetism turns 
into VBS order. In particular, the order of the transition is investigated using both 
finite-size scaling and the flowgram method of [l6|. 

4.1 Finite-Size Effects of {Wf) Near the Transition 

As we have seen in the previous section, antiferromagnetism is substantially weak- 
ened as the four-spin coupling Q increases. This manifests itself in the reduction 
of the staggered magnetization as well as in the increase of the characteristic 
length scale ^ = c/(27rps)- The higher order terms in the systematic expansion are 
suppressed as long as L 3> ^. In practice, this limits us to ^ ~ 10a which corre- 
sponds to Q/t/ ~ 5. As one approaches a second order phase transition, ^ diverges 
and the systematic effective theory is no longer applicable. Instead, in the vicinity 
of the phase transition, it is useful to employ finite-size scaling. 

In order to locate the transition it is natural to investigate the J / Q-dependence 
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of the spatial winding number squared (W^). In particular, in case of a second 
order phase transition, for sufficiently large volumes the various finite volume curves 
should all intersect at the critical coupling. Recently, such an analysis has been 



reported by Melko and Kaul [20j. We have verified explicitly that our Monte Carlo 
data are consistent with those of that study. In figure 5a we show a fit to those data 
for moderate volumes L/a = 32,40, and 48 near the transition using the finite-size 
scaling ansatz 

{W^) = f i^-j^L'"') =A + B^-j^L^'^ + O ■ (4.1) 

The fit is good, with x^/d.o.f. ~ 2, suggesting that the transition might actually be 
second order. In particular, the three finite volume curves intersect in one point, 
Jc/Q = 0.0375(5), and do not require an additive sub-leading correction CL~'^ 
to eq. fl4.1l) . This is consistent with Sandvik's earlier result obtained on smaller 
volumes L/a = 16, ...,32 which did require the inclusion of the sub-leading term. 
His fit led to ~ 2 which implies that the corrections are suppressed for large 
volumes. Remarkably, when Melko and Kaul's L/a = 64 data are included in 
the fit of eq. fl4.1l) . its quality degrades to x^/d-o.f. ~ 8. In fact, the L/a = 64 
curve does not pass through the intersection point of the smaller volume curves. 
Melko and Kaul attribute this behavior to sub-leading terms and finally come to 
the conclusion that there is a deconfined quantum critical point somewhere in the 
interval 0.038 < J/Q < 0.040. Indeed, a fit including an additional sub-leading term 
CL^^ is possible. However, the exponent uj is not well determined by the data. In 
order to obtain a stable fit, we have fixed uj to different values ranging from 0.01 
to 2.5, which all give more or less the same x^/d-o.f. ~ 2.5, but lead to different 



values of the critical coupling Jc. For example, fixing = 2, as suggested by [17 
one obtains Jc/Q = 0.0404(4) and z/ = 0.62(2). This fit is illustrated in figure 5b. 
On the other hand, when one fixes uo = 0.01, one obtains Jc/Q = 0.0438(7) and 
u = 0.62(2). Finally, when one excludes all but the largest volumes L/a = 40,48, 
and 64, a four-parameter fit becomes possible again. This fit, shown in figure 5c, is 
less good with x^/d.o.f. ~ 3.5 and it yields Jc/Q = 0.0398(6), which is inconsistent 
with the critical coupling obtained from the moderate volume data. Even larger 
volumes would be needed in order to decide if the curves will continue to intersect 
in the same point. 

To summarize, the moderate volume data (with L/a = 32,40, and 48) are well 
described by the four-parameter fit of eq. fl4.ip . while all data including those for 
L/a = 64 are not. These data can be described by a six-parameter fit including sub- 
leading corrections, but the data do not unambiguously determine the fit parameters. 
Since no sub-leading term is required to fit the moderate volume data, it seems 
strange that such corrections become necessary once larger volumes are included in 
the fit. We take this unusual behavior as a first indication that the transition may 
actually be weakly first order. 
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Figure 5: Three different fits of the spatial winding number squared (W^) as a func- 
tion of the coupling J/Q in the transition region. 
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J/Q = 0.038 



0.381 

20 30 40 50 60 70 

L/a 



Figure 6: Non-monotonic volume- dependence of {Wf) at JjQ = 0.038 near the 
critical coupling that may indicate a weak first order phase transition. 



Another observation that may cast some doubt on the second order nature of 
the transition is a non-monotonic behavior of (W^) near the transition, which is 
displayed in figure 6. Such behavior is typical for a first order phase transition. For 
example, in the VBS phase, at a point close to a first order transition, domains of an- 
tiferromagnetic phase can still exist. Thus, for small volumes, the antiferromagnetic 
domains may lead to a linear increase of (W^) with L. For larger volumes, the VBS 
phase will begin to dominate and (W^) will then decrease. This competition can 



lead to non-monotonic behavior. For these reasons, we think that the data of [17 
and jioj do not provide sufficiently convincing evidence that deconfined quantum 
criticality has actually been observed. Due to limited numerical resources, we have 
not been able to extend the analysis to substantially larger volumes. However, using 
a supercomputer this would definitely be possible and, in fact, highly desirable. In 
order to shed more light on the subtle issue of quantum criticality versus a weak 
first order transition, we now turn to an alternative method of analysis. 



4.2 Application of the Flowgram Method 

Kuklov, Prokof 'ev, Svistunov, and Troyer have developed a fiowgram method which 
is useful for distinguishing weakly first order from second order phase transitions 
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Iq . For our system, the flowgram method can be implemented as follows. We work 
on lattices of increasing size L at the inverse temperature given by (3Q = L/a. First, 
each individual spin configuration in the path integral is associated with either the 
antiferromagnetic or the VBS phase according to the following criterion. If all three 
winding numbers Wi, W2, and Wt are equal to zero, the configuration is associated 
with the VBS phase. On the other hand, if at least one of the three winding num- 
bers is non-zero, the configuration is associated with the antiferromagnetic phase. 
This criterion is natural, because in the infinite volume limit there is no winding in 
the VBS phase, while there is always some winding in the antiferromagnetic phase. 
One then defines a volume-dependent pseudo-critical coupling Jc{L) at which both 
competing phases have equal weight, i.e. the number of associated configurations 
is the same for both phases. It is important to note that, in the infinite volume 
limit, the pseudo-critical coupling Jc{L) approaches the true location of the phase 
transition both for a first and for a second order phase transition. The large vol- 
ume limit is now approached by simulating systems at the pseudo-critical coupling 
Jc{L) for increasing values of L. Defining the sum of the spatial and temporal wind- 
ing numbers squared as W"^ = + W2 + Wf, the quantity {W'^){Jc{L)) is then 
evaluated at the pseudo-critical coupling Jc{L). If the phase transition is second 
order, (iy^)(Jc(L)) will approach a constant for large L since ps then vanishes (i.e. 
^ = c/ {27rps) diverges) at the transition. On the other hand, if the transition is first 
order, with 50 percent probability the system still shows the characteristics of the 
antiferromagnet. Thus, for L ^ ^, {W^){Jc{L)) grows linearly with L. As we will 
see below, for 48a < L < 80a we indeed observe this behavior. 

We have implemented the Ferrenberg-Swendsen re-weighting method [23j in or- 
der to accurately locate the pseudo-critical coupling. Unlike in the rest of this paper, 
the simulations in this subsection have only been performed at two (instead of three) 
lattice spacings in discrete time. Both lattice spacings are close to the continuum 
limit and give consistent results. Instead of extrapolating to the continuum limit 
(which is less reliable with two than with three lattice spacings), in this subsection 
we quote our results at the smaller lattice spacing eQ = 0.05. A calculation closer to 
the continuum limit or, even better, using a continuous-time algorithm would still 
be useful. 

In order to investigate whether the phase transition is second or weakly first 
order, the values of {W'^) at the pseudo-critical coupling Jc{L) are illustrated in 
figure 7. For moderate volumes up to L = 48a the curve seems to level off, which 
would be characteristic of a second order phase transition. Indeed, as we have seen 
before, the moderate volume data for the spatial susceptibility are consistent with 
the finite-size scaling behavior of a second order phase transition. However, for larger 
volumes the curve rises linearly, thus indicating a weak first order phase transition. 
Of course, one cannot completely exclude that the curve may eventually level off at 
even larger volumes. We find this unlikely and conclude that our results cast serious 
doubt on the picture of deconfined quantum criticality painted in the earlier studies. 
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Figure 7: The sum of spatial and temporal winding numbers squared {W'^){Jc{L)) 
evaluated at the pseudo-critical coupling Jc{L) for increasing lattice size L. 



Lja 


24 


32 


48 


64 


80 


96 


Jc{L)/Q 


0.0311(4) 


0.0316(3) 


0.0337(4) 


0.0364(3) 


0.0384(3) 




J'c{L)/Q 


0.115(2) 


0.0871(4) 


0.0632(4) 


0.0544(5) 


0.0477(4) 


0.0445(4) 



Table 3: Values of the volume- dependent pseudo- critical couplings Jc{L) and J'^L) 
obtained with the Ferrenberg-Swendsen re-weighting method. 

Given the evidence for a weak first order transition, we like to determine the 
value of the critical coupling in the infinite volume. The values of the pseudo- 
critical coupling Jc{L) in a finite volume are summarized in table 3. Given the 
data for Jc{L) alone, it is non-trivial to extract the infinite volume critical coupling 
Jc = Jc{L — > oo). For this reason, we have defined another pseudo-critical coupling 
J'ciL), which also extrapolates to the correct limit, i.e. J^(L oo) = J^. In this 
case, we work at the inverse temperature given by (3Q = L/4a. Irrespective of the 
spatial winding numbers, if the temporal winding number Wt is equal to zero, the 
configuration is now associated with the VBS phase. On the other hand, if Wt is non- 
zero, the configuration is associated with the antiferromagnetic phase. As before, we 
define the volume-dependent pseudo-critical coupling J'^^L) such that both phases 
have equal weight. The values of J'^L) (again quoted at eQ = 0.05) are also listed in 
table 3. According to the finite-size scaling theory for first order phase transitions. 
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Figure 8: Fit of the pseudo-critical couplings Jc{L) (lower curve) and J'(.{L) (upper 
curve) shown as functions of a/L. 

using /9L^ oc L^, both finite- volume pseudo-critical couplings should approach their 
common infinite volume limit Jc as 

ML) = J. + C'^, .m = J. + C'^. (4,2) 

Interestingly, the two pseudo-critical couplings indeed converge to the same limit. 
A fit of Jc{L) and J'^L) to eg. (14.21) — shown in figure 8 — has x^/d.o.f. ~ 1.15 and 
yields the infinite- volume critical coupling Jc = 0.0396(6). Again, only the large- 
volume data show the characteristic behavior of a first order phase transition. It 
should be noted that the definition of J'^L) is less natural than the one of Jc{L), 
because it ignores the spatial winding numbers when configurations are associated 
with either of the two phases. In particular, J'^i^L) approaches the infinite- volume 
critical point Jc more slowly than Jc{L). Consequently, the ultimate large volume 
physics is more easily visible using the pseudo-critical coupling Jc{L). For example, 
the linear increase of {W'^){Jc{L)) with L, which sets in around L k, 50a, is not yet 
present in (VI^^)(J^(L)), and is expected to set in only on larger volumes. 
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5 Investigation of the VBS State 



As we have seen, the antiferromagnet is weakened and ultimately destroyed at a 
rather weak first order phase transition. Since the transition is so weak, at moderate 
volumes it is practically indistinguishable from a continuous transition. As a result, 
an approximate U{1) symmetry emerges dynamically as an enhancement of the 
discrete 90 degrees rotations of the square lattice. The other side of the phase 
transition is characterized by VBS order. However, the emergent U{1) symmetry 
makes it difficult to identify the nature of the VBS state as columnar or plaquette. 



5.1 Probability Distribution of the VBS Order Parameter 

In order to investigate the nature of the VBS order it is best to go away from the 
critical point as far as possible (assuming that no other phase transitions take place). 
In the following we thus work at Q/J = oo, which is obtained by putting J = 0. 
The corresponding probability distribution of the standard VBS order parameters 
p{Di, D2) has been determined by Sandvik for a 32^ lattice at zero temperature and 



it shows perfect U{1) symmetry [17|. The loop-cluster algorithm allows us to repeat 
this study for larger volumes, in this case using the probability distribution of the 
non-standard VBS order parameters p{Di,D2). As one sees from figure 9, even on 
a 96^ lattice at jSQ = 30 one does not see any deviation from the U{1) symmetry. 
Hence, our data do not allow us to identify the nature of the VBS order. 

At small Q/J, the loop cluster algorithm is extremely efficient with auto-cor- 
relations limited to at most a few sweeps. However, at larger values of Q/J, and 
especially at Q/J = 00 the algorithm suffers from a noticeable auto-correlation 
problem. This problem arises because the cluster algorithm, which is designed to 
update long-range spin correlations, can not efficiently shuffie spin-flip events from 
even to odd bonds. This causes slowing down in the VBS phase. Details concerning 
the algorithm and its performance will be discussed elsewhere. 



5.2 Antiferromagnetic Correlations in the VBS Phase 

In order to conflrm that antiferromagnet ism indeed disappears for large Q, we again 
consider Q/J = 00. We have simulated the staggered susceptibility as a function of 
the lattice size L. As one sees in flgure 10, at PQ = 50 the staggered susceptibility 
Xs increases with increasing space-time volume until it levels off around L ^ 50a. 
This shows that long (but not inflnite) range antiferromagnetic correlations survive 
even deep in the VBS phase. These data conflrm that antiferromagnetism is indeed 
destroyed in the VBS phase. However, again one needs to go to volumes larger than 
L ^ 50a in order to see the ultimate inflnite-volume behavior. 
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Figure 9: Probability distribution p{Di, D2) obtained on a 96^ lattice at (3Q = 30 and 
Q/J = 00. The observed U{1) rotation symmetry implies that we cannot identify 
the nature of the VBS phase as either columnar or plaquette. 
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Figure 10: The staggered susceptibility Xs in the VBS phase increases with increasing 
space-time volume until it levels off around L/a ^ 50 for j3Q = 50. 
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6 Conclusions 



We have employed a rather efficient loop-cluster algorithm to investigate the physics 
of the antiferromagnetic spin | Heisenberg model with an additional four-spin in- 
teraction. When the four-spin coupling is sufficiently strong, antiferromagnetism is 



destroyed and gives way to a VBS state. While Sandvik's pioneering study [17 



was 



limited to zero temperature and moderate volumes, just like the stochastic series 



expansion method of Melko and Kaul [20|], the cluster algorithm allows us to work 



at non-zero temperatures and large volumes. Using the cluster algorithm and ap- 



we 



plying the flowgram method of Kuklov, Prokof 'ev, Svistunov, and Troyer 
found numerical evidence for a weak ffist order phase transition, thus supporting the 
Ginzburg-Landau- Wilson paradigm. Interestingly, the same conclusion was reached 



in studies of the transition separating superfluidity from VBS order [IJ, |15|, 116 
The ffist order nature of the phase transition in the Heisenberg model with four- 
spin coupling Q implies that the idea of deconfined quantum criticality again lacks 
a physical system for which it is ffimly established. Hence, the proponents of this 
intriguing idea are challenged once more to suggest another microscopic system for 
which one expects this fascinating phenomenon to occur. 

It is interesting to ask why the phase transition separating antiferromagnetism 
from VBS order is so weakly first order. There must be a reason for the long 
correlation length, around 50a, even if it does not go to infinity. Perhaps, the ideas 
behind deconfined quantum criticality may still explain this behavior. 
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